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Previous dilepton measurements in heavy-ion collisions have mainly focused on invariant mass 
spectra to clarify in-medium changes of vector meson properties. However, a dimuon is characterized 
by two scales — invariant mass M and transverse momentum pr- Like transverse momentum 
spectra of hadrons, pr spectra of dileptons arise from an interplay between emission temperature 
and collective transverse flow, whereas the invariant mass is insensitive to flow. Having two control 
parameters of which only one is sensitive to flow allows at given M to characterize the emission 
region in terms of average temperature and flow. Thus, one is able to study what phases of the 
fireball evolution radiate into a given mass window. We demonstrate this technique using the dimuon 
transverse momentum spectra measured by the NA60 collaboration and present strong arguments 
that a thermalized evolution phase with T > 170 MeV leaves an imprint in the spectra. 
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I. INTRODUCTION 

The in-medium properties of vector mesons have long 
been considered a tool to observe chiral symmetry 
restoration experimentally. In order to provide the 
medium, decays of vector mesons produced in heavy-ion 
collisions have to be studied, and there is a long history 
of model calculations and measurements [if. So far, the 
experimentally accessible quantity has been the invariant 
mass M spectrum of dileptons, a quantity in which the 
in-medium spectral function, averaged over the evolution 
history of the medium, is reflected. 

With the NA60 high precision dimuon measurement [2] 
it has become clear that the averaged in-medium spec- 
tral function in the vector channel shows broadening of 
the p meson and no significant mass shift 0]. While it 
is difficult to establish this broadening as an unambigu- 
ous signal for chiral symmetry restoration, a suggestion 
has been made in that dilepton emission in the in- 
variant mass region above 0.9 GeV might be dominated 
by 477 processes if a considerable chiral mixing of vector 
and axial vector correlators is assumed. On the other 
hand, in we found that this region is dominated by 
partonic radiation. Thus, a tool to further characterize 
the emission region dominating a particular region in the 
invariant mass spectrum is clearly important before more 
definite conclusions can be drawn. 

Recently, NA60 has presented measurements of trans- 
verse momentum px spectra in three different mass bins 
0. This allows to study dimuon emission as a function 
of two scales M and pr with different properties. While 



A/ is invariant under Lorentz transformations and the 
invariant mass spectrum is hence (apart from detector 
acceptance effects) insensitive to collective flow of the 
medium, py-spectra are heavily influenced by flow. 
The importance of collective flow for the description of 
hadronic px spectra in heavy-ion collisions has been re- 
alized a long time ago (cf. e.g. [7|)- For experimental 
results from the Relativistic Heavy Ion Collider (RHIC) 
in Brookhaven, the analysis of both single hadron spectra 
and two-particle correlations in terms of thermal spectra 
subject to a flow fleld has been extremely fruitful to char- 
acterize the flnal state of the evolution Extension of 
the method to a parametrization of the full evolution in- 
stead of the flnal state has been demonstrated to give a 
good characterization of spectra and two-particle corre- 
lations at both SPS 158 AGeV flxed target Pb-Pb col- 
lisions [i| and RHIC 200 AGeV Au-Au colhsions 
In the following, we use the same techniques, adapted 
to the different nature of electromagnetic emission pro- 
cesses, to characterize the emission region dominating the 
PT-spectra in the different mass bins measured by NA60. 
Conceptually, this is similar to the ideas outlined in [11 1 
where elliptic flow is instead suggested to characterize 
the emission region. 



II. DILEPTON EMISSION AND FLOW 

Our model for dilepton emission in Inin collisions is out- 
lined in [sl , here we reproduce the essential details as far 
as needed to understand the flow analysis. 
We calculate the dilepton spectrum from the convolution 
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where the rate is found from the averaged virtual photon 
spectral function R{q,T, ps) as 



dN 



«2 R{q,T,pB) 



(2) 



Here, q represent the four-momentum of the emitted 
muon pair, T is the temperature of the emitting volume 
element and ps its baryon density. The finite Icpton 
mass rrLL can be accounted for by an additional factor 
of (1 + 2m\/lVp)^l - iml/A'P in Eq. Q. The fire- 
ball evolution encodes information on the evolution of 
radiating volume, temperature T and baryon chemical 
potential fiB, local transverse flow velocity vt, longitu- 
dinal rapidity 77 and in the late evolution stage chemical 
non-equilibrium properties such as pion chemical poten- 
tial /x,r and Kaon chemical potential pK- 
The model considers matter with a transverse entropy 
density s distribution 



s(r, r) = 1/ 1 + exp 



Rc{r) 



(3) 



where Rdr) = Rq + '-^t^ 
transverse flow field 



Onto this distribution, a 



Pt = atanh ?;j_(t) = r/Rci^r) ■ Pc{t) 



(4) 



is imposed with Pc{t) = atanh a^r. a±Tf = Vf is the 
final velocity at radius Rc{Tf) at freeze-out. Note that 
this is not the maximal velocity reached at the system's 
surface: The surface is dynamically calculated from the 
local entropy density and the equation of state via the 
condition T = Tp- At r/, a thermalized system exists 
only close to i? = and Rc{Tf) is outside the thermalized 
region. Since the position of the Cooper-Frye surface 
changes in a non-trivial way the model cannot completely 
be characterized by one number for the flow. We have 
also investigated flow fields with dependences p ^ and 
P ^ which are suitably normalized such as to lead 
to hadronic pr spectra with the same slope as the form 
Eq. (dl (see [3, [101 for details about the calculation of 
hadronic emission from the model). 
Let us recall the role of transverse flow: The factor p^w^ 
in Eq. Q evaluates for a process at midrapidity to 



E ~ niT cosh(pT) — Pt sinh(pT) cos((/)) 



(5) 



where <j) is the angle between flow direction and dilcpton 
momentum and E the energy of a virtual photon mea- 
sured with M,pt, rriT ~ \/ A/^ + p\. as it appears in the 
local restframe of a fluid element with radial rapidity p^- 
For a particular choice of flow field and transverse density 
distribution and a Boltzmann distribution instead of the 



rate 



R[q,T,pB) 



Y in Eq. ^ (as appropriate for hadrons), 
the spatial integral of the convolution of the flow field 
with density yields analytical expressions for the spec- 
trum 0, Q , however we solve the integral numerically. 



It can be seen as a general property of Eq. ([5]) that the 
resulting effect on the spectral slope will in general be pr 
dependent. This is apparent from the two limiting cases 
Pt ^ M where the equation becomes E k, M cosh p and 
Pt 3> M where one finds E « p7'exp[p] for (jj = n and 
E Pt exp[— p] for = 0. It is well known that the effect 
of fiow makes the spectral slope harder, in the case of 
large masses and flow values even leading to a flattening 
of the spectrum, the 'knee' seen e.g. for protons Q. This 
hardening of the slope is mass dependent as long as pT 
is of the order of M . 

The slope of a spectum in the model thus arises as follows: 
The temperature T in Eq. ^ sets the basic scale possibly 
modified by the momentum (and density) dependence of 
the spectral function R{q,T, pb)- We point out that the 
modifications of the p^-slope due to this dependence are 
small for all investigated thermal sources, namely for the 
in-medium hadronic sources (in-medium vector mesons 
and four-pion annihilation) as well as for the partonic 
source. 

In the absence of fiow, a dilepton pT spectrum at given M 
arises from a superposition of contributions from different 
T, weighted by the volume radiating into the acceptance 
at this temperature. Since T drops as a result of volume 
expansion, this means that the low T spectra will typi- 
cally contribute more to the pt integrated yield as long as 
M is sufficiently small, thus adding all contributions the 
low Pt region will be dominated by late, low T emission 
whereas at larger pT the harder slopes of earlier emission 
from hotter regions will gradually become visible. 

The presence of a flow field changes this picture some- 
what: Since T and flow are anticorrelated in the model 
(initially T is large and pT small, at late times the rela- 
tion is reversed), there is no cold emission source with- 
out flow distortion in the model, but hot emission re- 
gions are unaffected by flow. For a single temperature, 
the flow effect leads to a pT-dependent hardening of the 
spectrum. This blueshift becomes maximal at some pt 
(dependent on the values of M and pt), the 'knee', and 
then drops again. The final pT spectrum arises as the 
integral over all volume elements with different temper- 
atures and blueshift values, i.e. the knee position is 
averaged and the pT-depcndent sequence of hardening- 
softcning of the slope indicating a 'knee' becomes invis- 
ible. The general feature that the soft slopes dominate 
at low Pt whereas the hardest slope becomes visible at 
high Pt returns instead (this is, in fact, the major dif- 
ference of e.m. spectral slopes as compared to hadronic 
spectral slopes where only contributions from a narrow 
temperature region need to be considered). 

The model also imposes a longitudinal flow field, however 
since it has little relevance for the slope of pt spectra we 
do not discuss it further, the interested reader is referred 
to 1]. 
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III. COMPARISON WITH NA60 DATA 

The NA60 collaboration has measured px spectra in 
three different invariant mass windows: A low mass win- 
dow with 0.4 GeV < M < 0.6 GeV, the p-hke region 
with 0.6 GeV < M < 0.9 GeV and the high mass region 
1.0 GeV < M < 1.4 GeV. While our model as outhned 
in has predicted the general systcmatics and order 
of magnitude of the pT slopes surprisingly well 0] just 
based on scaling relations from Pb-Pb collisions to the 
smaller In-In system [Bj, in the following we do the fol- 
lowing minor readjustments to improve agreement with 
data: For the flow field, we use a dependence p ~ ^/r 
instead of Eq. (j4|) . In addition, we increase the (unphysi- 
cal) parameter Vf to 0.6c such that the maximum flow at 
reached at the medium surface close to r/ is about 0.57 
c. The overall change in the spectral slopes as compared 
to the results presented in [1,0] are less than 5%. 
Let us start the discussion with the p-like mass region. 
Here, the invariant mass spectra show a clear peak with 
vacuum width which has been identified with p mesons 
decaying after thermal decoupling of the hot matter in 
vacuum 0, [l^l- If so, the spectral slope of this contribu- 
tion should coincide with the spectral slope of a hadronic 
spectrum at the p mass. However, in both analyses (sl.lT^ 
the vacuum contribution does not dominate the yield in 
the invariant mass window. We may however note that 
the vacuum contribution is expected to reflect a maximal 
amount of flow (and hence hardest slope), as it arises 
from the three-dimensional boundary region of the ther- 
malized four-volume and the flow field is always largest 
at the edges. Thus, even if this contribution does not 
dominate when integrated over px, at sufficiently high 
Pt it should. 

If so, the slope at sufficiently high px can be estimated 
using a phenomcnological formula from [isj . There it 
was shown that the systematics of hadronic spectra at 
different mass can be described roughly by a fit of the 
form 



1 dN 
rriT drriT 



A exp 



rriT 



(6) 



and that the pT spectrum in this region approaches the 
expected value with increasing pT as a strong indication 
that the data indeed show the vacuum p. A decompo- 
sition into different sources of the theoretical prediction 
of the pT-spectrum integrated over 0.6 GeV < M < 0.9 
GeV mass region indeed shows that the vacuum-p con- 
tribution is the dominant source for pt ^ 1 GeV. This 
puts severe limits on the flow field close to r^, and by 
backward extrapolation of the transverse and longitudi- 
nal flow fields also on somewhat earlier times. 
We present the full calculation of both spectrum and ef- 
fective temperature T* as extracted from a local fit of the 
spectra to the r.h.s. of Eq. ^ in all three mass windows 
compared with the data in Fig. [TJ As can be seen, the 
calculation with the flow field constrained by the vacuum 
rho decay describes the spectra in general well except in 
the low Pt region with pT < 0.5 GeV. 
Let us discuss the three mass windows in turn: Both the 
low mass and the p like region receive their dominant 
contribution in our model from the late hadronic evolu- 
tion phase, i.e. from regions with temperatures T < 170 
MeV. This implies that the blueshift of the spectra by 
flow is large — it amounts to about 70-80 MeV and is 
hence even substantially larger than the variation of the 
actual temperature in the hadronic phase which is only 
40 MeV. In particular, it is substantially larger than the 
uncertainty introduced by the choice of the proper px de- 
pendence of the spectral function in the hadronic phase 
— even a drastic change like neglecting any pT depen- 
dence of the spectral function does not lead to more than 
10 MeV change in slope. 

The difference between the low mass and the p like-region 
is essentially caused by the contribution of the vacuum p 
which receives the maximum amount of flow (see above) . 
In contrast (as also apparent from T*) the high mass bin 
is in our calculation dominated by partonic radiation. 
This predominantly takes place before the transverse flow 
fleld grows large and hence T* in this momentum win- 
dow approximately reflects the true temperature of the 
system. The model results for T* as a function of M are 
shown in Fig. O The prominent role of the vacuum rho 
contribution is clearly apparent from this figure. 



where 



IV. THE LOW PT REGION 



T*=T + m,{vTf (7) 

with rrii the mass of the hadron in question, T* the appar- 
ent spectral slope and T the true temperature of the emis- 
sion region at freeze-out. Within reasonable limits, this 
formula reproduces the more microscopic considerations 
of e.g. Ji] quite well. Averaging our flow proflle over the 
fireball breakup time and evaluating this relation, we find 
T* « 0.265 MeV, i.e. a rather hard spectrum. However, 
this is in fact not far from the experimentally observed 
value, given the simple estimate we have made. We take 
the fact that a peak with vacuum p width is observed 



As apparent from Fig. [U the model fails to reproduce 
the low Pt region in all three mass bins. Experimentally, 
values of T* as low as 100 MeV are observed. Recall- 
ing that applying a flow fleld to a pt spectrum leads to a 
blueshift, i.e. higher apparent temperature, we may thus 
characterize the source radiating into the low pt region 
as having a temperature T < 100 MeV and no strong 
flow. The strict anticorrelation of T and pT in any hy- 
drodynamically motivated model makes it appear very 
questionable that this could be a thermal source. A sim- 
ilar softening of the low pT slope is observed e.g. in pion 
spectra and could there be ascribed to resonance decays 
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FIG. 1: Left panel: pr spectra in three different mass windows as obtained by NA60 @| as compared to the model calculation. 
Right panel; Effective temperature of the spectra determined from an exponential fit for a pr window of 0.8 GeV (0.6 GeV for 
the first data point) as a function of the window center. 
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FIG. 2: Spectral slope T* at pr = 1.5 GeV calculated as a 
function of invariant mass M. 

, although it is at present unclear what decay process 
could be responsible for the source dominating the low 
Pt region in the dimuon spectra. 

Since it can clearly be demonstrated that a thermal de- 
scription of a flowing medium cannot possible describe 
the source properly, we will not address the region of 
Pt ^ 0.5 GeV in the following discussion. 

V. THE LARGE M REGION 

While the data seem to support our idea that the mass 
region above 1 GeV is dominated by partonic radiation, 
the question remains if a hadronic description would not 
likewise be possible. We will address this question by 
requiring that the correlation of px and T must be kept 
such that the low mass and p-like region are described 
properly. Given this constraint, we ask what source char- 
acteristics are compatible with the measured slope. 
Note that the fact that T* is lowest in the high mass re- 
gion is far from trivial: If the source radiating into high 
masses were characterized by similar average tempera- 



ture and flow as the source in the other mass bins, the 
slope would be expected to be hardest as the blueshift of 
spectra in a flow field increases with mass (see Eqs. (|5I7|) ). 
In fact, if the source had the same average emission tem- 
perature, we would expect T* = 250 — 260 MeV given 
our flow field. Thus, it is evident that there must be a 
transition to a qualitatively different source. 
Assuming that the dominant source is hadronic, T < 170 
MeV, but not close to Tp (see above) we can next ask 
how big a window in T with associated changes in the 
flow field one can allow before the low mass and the p-like 
mass bins cannot be described any more by the model. 
We find that the region from 170 down to 160 MeV is in 
principle compatible with the data, the region down to 
150 MeV is marginally compatible, but the source cannot 
dominate the high mass yield for T < 150 MeV, otherwise 
the strong flow field (required by the lower mass bins 
which are dominated by emission from this temperature) 
pushes T* above the data. 

Thus, if four-pion annihilation processes were the domi- 
nantc source like suggested in [3l they would have to cease 
contributing rather suddenly around T ~ 160 MeV. 

Employing this argument solely based on the p^-slope of 
the di-muon spectrum, one cannot rule out that the high 
mass radiation comes dominantly from a hadronic source 
with temperatures in the range of 150 MeV < T < 170 
MeV, although a substantial amount of fine-tuning of the 
source is required for this assumption. However, assum- 
ing that such a fine-tuning is performed, the measured 
slope cannot be described by a hadronic source without 
flow, instead a substantial flow of order {vt) > 0.2 is 
needed. Such transverse flow indicates a precursor state 
which cannot be a mixed phase of a 1st order phase tran- 
sition (as no flow is built up from a mixed phase). Thus, 
even if the radiation in the high mass window would be of 
hadronic origin, its slope still points towards a partonic 
evolution before. 

Further constraints can be derived if the fact is taken 
into account that four-pion annihilation can have strong 
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contributions to the yield from emissions at temperatures 
close to Tf = ISOMeV < 160 MeV due to the growing 
fugacity factor cxp[A^Tr{T) /T] with smaller T. In an ex- 
plicit calculation with such a dimuon source, we employed 
thermal four-pion annihiliation rates as calculated in the 
approach This approach is constrained e.g. by the 
most recent BaBar data [l5| on the inverse process in 
the vacuum, namely lepton annihilation into four-pions. 
The calculation indeed shows that even if fine-tuning of 
the source is performed to enforce compatibility of T* 
with the pT'-slope data in the higher mass window, the 
di-muon yield from these processes is not enough to ac- 
count for the absolute yield as observed in the mass re- 
gion for M >1.0 GeV assuming these were the dominant 
source p^ . 

We point out that this investigation only implies that the 
dominant source in the M >1.0 GeV region cannot be 
four-pion annihilation processes. It is still possible that 
these processes contribute subdominantly to the M > 1 
GeV dimuon yield in that region of the mass spectrum. 
A careful analysis of the quantitative strength of such 
a possible subdominant contribution of four-pion anni- 
hilation processes to the spectrum has also to account 
for the fact that multiparticle annihilation processes can 
be expected (due to their dependence on the pion 
density p^) to decouple earlier than two particle anni- 
hilation. Even taking fugacity factors into account, 
in the calculation drops by more than a factor two from 
Tc to Tp. While common modelling assumptions treat 
all processes in kinetic equilibrium down to the same de- 
coupling temperature Tp, the experimental findings of 
momentum equilibration can be accounted for by 2 2 
processes alone, and more microscopic descriptions of the 
freeze-out process (e.g. hadronic cascades as in show 
a sequencial decoupling of different processes. Assuming 
full equilibrium for multipion annihilation down to Tp 
thus constitutes an upper limit for the expected contri- 
bution. Further detailed analysis of multipion annihila- 
tion is beyond the scope of the present paper, and will 
be presented elsewhere [l6j . 

VI. SUMMARY 

We have performed an analysis of the pt spectra by 
dimuons based on a picture of thermalized matter ex- 



panding with a given flow field that has been rather suc- 
cesful in the analysis of hadronic spectra. The essentail 
features of the spectra can be understood from thermal 
scale and blueshift by flow, thus our results do not cru- 
cially depend on details how hadronic or partonic emis- 
sion is modelled. Requiring that the same flow field and 
the same pT — T anticorrelation underlies all relevant 
sources, we have argued that the high mass region 1.0 
GeV < M < 1.4 GeV is most likely dominated by a 
partonic phase. 

While considerable fine-tuning of the model in princi- 
ple allows for a hadronic source at T close to the phase 
boundary Tc in this region, the implication is still that 
substantial flow is needed to blueshift such a source close 
to the data which can only arise from a pre-hadronic 
phase different from a mixed phase. In a thermal- 
ized description, this implies partonic evolution. We 
furthermore have explicilty excluded four-pion annihila- 
tion contributions as such a dominant hadronic source 
since its dimuon yield would not be sufRcicnt to account 
solely for the observed yield in the mass spectrum above 
M > 1 GeV even if the fine-tuning to account for the pr- 
spcctra's slope is performed. This does not exclude four- 
pion annihilation as a subdominant component in the 
spectrum. An interpretation of the high mass dilepton 
radiation in terms of a dominant partonic source seems 
more natural, as a partonic evolution before hadronic ra- 
diation is required by the data in any case. 
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